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INTRODUCTION 

The rapid advances made in the field of autovatic controls since 
Worid War II have found more and more applications in the area of ship 
propulsion. One of the principal. driving forces behind this is the 
desire to operata plants safely and at or near thoir "optimal" operating 
point. The safe coordinated operation of some ship power rlants cone 
sisting of hizh performance prime movers, sophisticated couplings, and 
propellers is no trivial matter. This is especially true during periods 
when fast ship response is desised, o.ge, emergency stopping or acceler= — 
ating. 

To be able to control a propulsion plant, on2 must bs able to 
model it and knew how it behaves under various disturbances or changes 
in inputs. The alternative is a costly experimental program for each 
system. The problem for the control engineer now becomes one of developing 
or choosing a model that suitably describes the characteristics of the 
particular power train couponents. These component models may be either 
physical or functional. A physical model may accurately describe a 
component's characteristics using, Zor example, Newson's Laws of Motion, 
thormodynamic relationships, conservation of mass relations, etc. These 
provide a valuable insight of the physical precesses occurring in the 
operation of a component, and enatle one to perceive what simplifications 
can be made to a model. Some components, on the other hand, may prove 
to be very non-linear and to describe their operational characteristics 
using physical laws may result in pabneda nies expressions which are not 
suitable for use in a simulation study. In this situation a simpler 
mathematical expression is developed which closely approximates the 


component's characteristics over a reasonable range of inputs. 





INTRODUCTION (Cont'd 

The philosophy behind the study described in this thesis has been 
to keep the models simple, but yet reasonable. To quote P. J. Kiviat, 
who prepared a RAND memorandum on digital simulation modeling concepts, 

"the model should only be as detailed as is necessary to 

answer the question at hand; it should be so designed, 

however, so it can be expanded to include more detail 

without inordinate cost in those model areas which have 

a high probability of becoming subsequent subjects of 

concern.''(1) 

Once the component models are selected, drive train models can be 
synthesized and prepared for simulatione No discussion will be made 
concerning the relative merits of digital simulation compared to analog 
or analogedigital (hybrid) simulation. The simulation literature is 
replete with these discussions. Digital simulation was decided at the 
outsed. 

A major part of this study has been concentrated on writing computer 
programs for prime movers, propellers, and transient operations. These 
programs were cevelcped to investigate various control schemes for marine 
power plants. 

The gas turbines selected for simulation are the Pratt and Whitney 
FTA~2 and the General Electric LM2500-A. <A medium speed and a low speed 
diesel were also selected. The medium speed diesel selected is the 
Mirrless KV Major 12; the low speed diesel selected is the Burmeister 


and Wain 7K98FF. The main reasons for selecting these particular engines 


were the availability of suitable data at reasonable power levels. 





1. 


GAS TURBINE MODELS 


lLel Introduction 


Gas turbines, as applied to ship propulsion, usually connote 
the inclusion of a power turbine or "free" turbine, which 
provides the driving torque. The power turbine is driven by 
the hot exhaust gases of the "gas generator." Thus when 
marine engineers refer to gas turbines, they usually include 
the combination of gas generator and power turbine in their 
meaning. 

Tnere are, of course, exceptions to the statement cone 
cerning separate power turbines. The Koelin Class escort 
frigates built for the German Federal Navy are driven by 
gas turbines which do not have a separate power turodine. 

{ne gas generator is coupled through - clutch to the drive 
traine(2) 
Aireraft Gas Turbine Models 

- During the literature survey part of this study, several 
gas turbine models useful for aircraft propulsion simulation 
were found. In Reference (3), for example, Cottington pre- 
sents a simulation study of a two-spool turbojet engine, and 
also covers the design of a digital three-term controller 
using digital simulation methods. [In trying to determine an 
"optimum" set of gains, he treated the high pressure rotor 
system as a critically damped second order system. Values 
of these optimized gain settings for varying values of high 
pressure rotor speeds are presented. 

Mueller has developed a single spool gas turbine model 
from one-dimensional flow relationships.(4) He assumes that 


the only important dynamics are those of the rotating compressor- 





le GAS TURBINE MODETS (Cont'd) 
1.2 Aircraft Gas Turbine Models (Cont'd) 
turbine spool. Mueller justifies this by noting that gas 
disturbances travel through the engine at the speed of sound, 
whereas a spool acceleration or deceleration takes some seconds. 
No results of his simulations are presented. 

In References (5) and (6), models are derived from a 
linear analysis of gas turbine dynamics. The model presented 
in the latter reference is highly detailed and is developed 
Using nonedimensionalized parameters. To use this model, 
one must have access to the machinery performance maps. The 
prominent feature of this model is that it assumes stepwise 
Linearity. The compressor performance map is utilized to 
obtain the coefficients in the linearized equation at each 
operating point. 

In References (7) and (8), Saravanamuttoo and Fawke dise 
cuss the application of both analog and digital computer methods 
of simulating gas turbine dynamic performance. They approached 
the problem from the viewpoint of engineering thermodynamics, 
using the normal compressor and turbine characteristics. 

Figure (1-1) shows the results obtained by Saravanamuttoo and 
Fawke for a single-spool turbojet simulation, both digital 

and analog, along with the results of actual engine tests. 

It can be seen that the acceleration predicted by both simue 
lations was much faster than the actual engine. An examination 
of the test results on the engine showed that there was a lag 


in the rise of turbine entry temperature following a fuel change. 





( RPM X* 1000 ) 


SPEED 


Tae TTT 8 
pen Gn Iigieey 


ENGIN 


we Ce iy ie 
O30 es V 


FEVe Vet rey 
Lar at 


bah} 
~, 
wa 


TIME (SECONDS) 


Seer ef V SaNsleeetca Gs re | SD in oir he 
lie LCO--! Cs ok oer uss LUW YD CC 


Fige (1-1) 


4A 


ce eae 


LD 





- oF reser 


’ 
ta Sack 


a 





1. 


GAS TURBINE MODELS (Cont'd) 


1.2 Aircraft Gas Turbine Models 





Saravanamuttoo and Fawke's model, however, had assumed 
immediate response, although they realized there would be 
some lag due to the dynamics of the fuel injection system. 
They decided to simulate this lag by postulating that after 
@ sudden change in fuel flow, the combustion efficiency drops 
instantaneously and then recovers exponentially. This drop 
in efficiency was made proportional to the percentage increase 
in fuel flow. The constant of proportionality was determined 
such that the turbine entry temperature before and after the 
transient was unchanged from the steady state values. The 
time constant for tha recovery of conbustion cfficiency was 
taken to be inversely proportional to the air flow through 


the combustor, ieGe; 
a , Wo 
Z= C, W 


T - time constant for recovery of conbustion efficiency 

Z,- time constant at design speed 

W,- air flow at design speed 

We air flow 

To achieve flow compatability in their digital model, 
Saravanamuttoo and Fawke added lags in the propagation of 
pressure signals through the engine. The lags were introe 
duced by considering the engine components to be separated 
by specified volumes; these volumes produced the desired 


capacitive effect, 
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1.4 


GAS TURBINE MODELS (Cont'd 


Marine Gas Turbine Models 

During the early parts of the literature survey, a Russian 
title on analog and digital models of marine gas turbines 
was found in a gas turbine newsletter. This book, listed 
as Reference (9), was written by Bichayev and published by 
the Shipbuilding Press of Leningrad in 1969. Although it 
is still untranslated, the book appears to be a major work 
in the field of marine gas turbine simulation, judging from 
the many analog computer diagrams, compressor maps, and 
equations. The cost of having it translated (over $900) 
was prohibitive. 

An unedited rough translation of a Russian work by Lurye 
discusses mathematical modeling of a marine gas turbine engine 
with a real fuel apparatus.(10) lLurye's model is presented 
in the form of linearized differential equations. The equations 
describing Lurye's model are reproduced in Appendix A. 

Using Saravenamuttoo's technique, one can extend his 
model to a gas turbine having a free turodine. This is also 
presented in Appendix A. 

Torque.Speed. Characteristics Model 

Ce Je Rubis has developed a method based on the free 
turbine torque-speed characteristics to model a gas turbine's 
action.(11,12) The free turbine torque-speed characteristics 
can be derived from the steady-state power.-speed characteristics 
which show lines of constant SFC. The corresponding curves 


for the Pratt and Whitney FTA-2 and General Electric LM2500-A 
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1.4 Torque~-Speed Characteristics Model (Cont'd 


GAS TURBINE MODSLS (Cont'd) 





are given in Appendix A. The torague-speed characteristics 
have been used here as the model selected for simulating the 
performance of the marine gas turbine. 

The unrefined model, used to describe the transient 
performance of the turbine operates in the following manner. 
Referring to Figure (12), assume that the initial steady-state 
operating point is at A operating at torque To and speed Nos 
and assume that the engine is to accelerate to point Be In 
a time interval At, fuel is increased by AF. This increase 
in fuel produces an instantaneous change in torque AT such 
that the new torque is Ty + AT. This increase in torgue will 
result in a change in power turbine rotation AN along the 
constant fuel rate line Fo +4QF. The new starting point for 
the next time interval is at the intersection of the R,+AF 
constant fuel rate line and the new rotational speed No +AN. 


The dotted lines illustrate the paths taken during the time 


steps, while the heavy line is the torque path. 


1.5 An Improved Model 


The methodology of this model is not new. The torque- 
speed characteristics have been used in the past to model 
diesel engines. The idea of obtaining instantaneous torque 
for a change in fuel rate is much closer to the real thing in 
a diesel than it is in a gas turbine. Comparisons with engine 
tests have shown that the predicted values of torques cbtained 


from the gas turbine characteristics are higher than the 
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1. GAS TURBINE MODELS (Cont'd 


1-5 An Improved Model (Cont'd) 


measured values. The differences in torque and the resulting 





differences in speed were measured by Rubis who conducted 
engine tests with the FTM4A-2. He also used the torque-speed 
curves in a simulation between the same two power levels. 
Figure (1-3) shows how the simulation results in a higher 
acceleration because of the higher torques. This graph should 
be compared with figure (I-/). 

It is evident from the comparison of figures (i-/) and (3) 
that what is needed to improve this mode] is a correction 
factor, which could have a physical explanation similar to 
that offered by Saravanamuttoo and Fawke for their model; 
namely a reduction in combustion efficiency. The results of 
other transient tests conducted with the FTA-2 were presented 
by Rubis in Reference (13). In this paper, Rubis presents 
correction factors to the torque-speed curve model for engine 
accelerations. This curve is reproduced herein as figure (1-4). 
One can observe that the curve for accelerations from idle can 
be expressed as a decaying exponential similar to the proposal 
by Saravanamuttoo and Fawke, vize, 

Ke = .3 + 7 aes ) 

This expression, which creates a form of delay or lag in the 
transient torque produced by the model, was included in the 
simulations in this study. It should be clear that this lag 
is applied to the transient toraue obtained from the model 


and not to the whole torque. 
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1. GAS TURBINE MODELS (Cont'd) 


1.5 An Improved Model (Cont'd) 


The model appears to be acceptable for ship propulsion 
studies because it is simple; it also permits loading other 
characteristics of the ship drive train into a medium sized 


computer without exceeding its core size. 





2. DIESEL ENGINE 
2el Introduction 


Lak 


Tne traditional use of steam power plants to propel 
large American built ships has created comewhet of a void 
in this country's technical literature on diesel ship pro- 
pulsion. The successful use of diesel engines for ship 
propulsion by foreign shipbuilders, on the other hand, has 
been an established fact long before the first digital com- 
puter was ever built, and long before analog and digital 
Simulations became established as accurate techniques for 
modeling and testing ship propulsion plants prior to their 
construction. Consequently there has been very little demand 
from industry, either foreign or domestic, for diesel plant 
simulations. However, some published data were obtained 
and used in this study. Since the use of high speed diesels 
has been fairly restricted to small craft, the simulations 


have been limited to medium and low speed diesels, 


Presentation of Models 


Woodward presents in Reference (14) a simple, straight- 
forward treatment of the diesel engine as applied to ship 
propulsion. A widely used expression given in Reference (14) 
and other engine textbooks, is the formula for indicated 


horsepower (IHP): 


_ PLAN 
THP = 330650 (HP) 


P = indicated mean effective pressure (1bs/in*) 


L = stroke (ft) 
re eee " ae ; 
A = total piston area (in”) for & mylti-cylinder engine 
: RPM + ; 
N = RPM (twoestroke engine), at (four-stroke engine) 
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Presentation of Models (Cont'd 

This power is that produced in the engine cylinders, and is 
greater than that available at the output shaft because of 
friction and other internal loads (e.g, water pumps, fuel pumps, 
turbochangers, etc.). What is available at the output shaft 

is termed brakehorsepower (BHP). The ratio of brakehorsepower 
(BHP) to indicated horsepower (IHP) is the mechanical efficiency 


y) of the engine, i.e., 


. DEP 
Rm ~ TEP 


or BHP = 7 PLAN 
33000 

The product Nae is called the brake mean effective pressure 
(BMEP). The indicated mean effective pressure, P, can be 
directly measured and is proportional to the energy released 
by combustion in the cylinder. | This energy is proportional 
to the amount of fuel injected in each power stroke. Most 
diesels utilize positive displacement pumps for fuel injection. 
The amount injected per power stroke for this type of pumping 
system should be the same. The stroke, L, and the total piston 
area, A, are, of course, fixed for a particular diesel. If 
one assumes constant mechanical efficiency, then it can be 
seen that for a particular fuel setting, BHP is a linear 
function of N, 1ieCe; 

PHP = KN 


K = NA for a particular fuel injection setting 
33900 
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2. DIESEL ENGINE (Cont'd) 


2Ze2 Presentation of Models 





Since rotational power is the product of torque, ©, and speed, 
t&, it then follows that for a two stroke engine, 

cy= 20 N (rad/min) 

Tw= C27 N  (ft-lb/min) 


oe 


i K as previously defined 

20 

i1.€., toraue is constant for particular fuel injection setting. 
The above is known as the ideal torque-rpm and povwer-rpm 
diesel relationships, and are graphically illustrated in 
Figure 2el and 2-2. 

The important assumption of this model is the one concern-~ 
ing constant mechanical efficiency. How close does this 
approach the actual case? In Reference (15), Pounder devotes 
an entire chapter to the Burmeister and Wain Engines and states 
a mechanical efficiency of 91% in his discussion on B &W low 
speed diesels. (Note: this figure was confirmed by a representa- 
tive of the Burmsister & Wain Corporation on a recent visit to 
MIT). 

Jt was indicated before that the indicated mean effective 
pressure is proportional to the fuel injected rate. An exe 
tension of the ideal model then is to say that torque is not 
only constant for a given fuel setting, but also proportional 
to that fuel setting. Jn other words, torque is a linear 
function of fuel setting. Obviously, other extensions can 


be made such as relating toraue to fuel rack position or to 





1a 


Ze DIESEL ENGINE (Cont'd 
202 Presentation of Models (Cont'd) 
the novement of some fuel control actuator. 
Lewis, Lecourt, and Scoville used a variation of tho 
ideal diesel engine in Reference (16). The torque developed 
within the engine was designated One and the friction torque 


Qre? The output torque is the difference between the two 


e 


torques. was postulated as linearly related to fuel 


Spe 
rack position, 2, 1-Ge, 

Qn = ‘Ke 
and the friction torque was taken as linearly related to engine 
speed, 

Qre = ho +A, Ne 
It can be seen that this model does not assume constant 
mechanical efficiency. 

In Reference (17), Garvey presents a linearized differene- 
tial equation describing the motion of a diesel engine. He 
assumed the output torque developed by the engine to be pro- 
portional to governor position, 1-€e, 

T = Cz 
T « chanze in torque 


Z governor position 


c « governor gain 


Parker uses the same model with a delay (18), 


T = ae 


ab Toraue 


c¢ = governor gain 


N 
il 


governor position 





13 


2. DIESEL ENGINE (Cont'd) 


2e2 Presentation of Models (Cont'd) 


Z ott 





s = differential operator 

t = time delay 
Low Speed Diesel Model 

The ideal diesel engine model in conjunction with a delay 
in the form of an exponential lag was chosen as the method 
to simulate a Burmeister & Wain 7K98FF low speed diesel. 
In Reference (19). Meuwlengracht presents the results from 
test-bed trials of a 7K98FF. From these results, and assuming 
the ideal diesel characteristics, one can get a brakehorse- 
power-speed curve and from this curve, derive the flat torque 
speed cheracteristics for constant fuel settings. These 
curves are all included in Appendix B. 


Mediun Speed Diesel Model 









Reference (20) describos the development of high per. 
formance medium speed diesel, the Mirrless KV Major 12. 
Since the fuel map for this diesel is included in this article, 
it was selected as the diesel to modele From the fuel map, 
and using the elementary relaticns presented earlier, one 
ean derive the engine's torque spcea charucteristics for cone 
stant fuel settings. In the process of drawing these charac. 
teristics, it was observed that torques could be approximated 
quite well as being proporticnal to fuel setting. Thus 
the expression, 


i 
x 


Kx (torque) 
fuel setting 


HH) 


i) 


K = proportionality constant 


was used in the simulations. 





oe 


LY 


DIESEL ENGINE (Cont'd) 


2-5 Delays 


26 


As in the gas turbine, the diesel engine does not 
instantaneously produce a change in torque with an increase 
in fuel ratae For a four-stroke cycle multi-cylinder diesel 
engine, Parker recommends a valve of =2 when N equals engine 
rp; this is an epproximation for diesel engines with two 
or more cylinders firing per revolution.(21) For the more 


general case, Reference (16) recommends 


2 < Thmy<B + 2 
Re Re ReP 


Nee? designed engine speed 

Pp : number of cylinders firing per revolution 
Discussion 

The torque-fuel setting characteristics were found to 
be linear for the medium speed-diesel but not ouite linear 
for the low speed engine. The characteristics for the low. 
speed diesel were therefore stored in the computer. For each 
diesel, a change in fnel setting results in @ change in torque. 
This torque was tempered by a delay in the form at a decaying 
exponential. This method is much simpler than that employed 
for the gas turbine models. since diesel torque is independent 
ox speed in the enginets operating range. 

For ship propulsion studies being considered here, it was 
not considered necessary to seek a model which would produce 
the oscillatory torque that is characteristic of internal ccom- 


bustion engines (cf. Reference 16). 





20 


3e CLUTCHES AND HYDRODYNAMIC TORQUE CONVERTERS 


3.1 Introduction 


The selection of a fixed pitch propeller and a none 
reversing prime mover will necessarily mean the incorporation 
of a clutch or torque converter in the drive train. Gas 
turbines, for example, are limited to undirectional movement; 
and some diesels, on the other hand, are designed for operation 
in two directions. It will be observed later that clutch 
engagement and disengagement times play a critical part in 
the emergency stopping of a ship. The selection of a proper 
clutch for a drive train thus is no trivial matter. In the 
same vein, the selection of a suitable clutch model to simu- 
late the emergency stopping of a ship is important in order 
to avoid erroneous results. 

Fajk Marine Airflex Clutch 

Figure 3-1 shows the principal part of the Airflex clutch. 
The clutch is aeeical ty operated drumetype dry friction 
clutch consisting of an annular rubber/fabric air gland bonded 
to an outer, flanged steel rim and lined with friction shoes 
which engage an inner iron drum when the gland is inflated. 

A never air-cooled version utilizes extruded aluminum shoe 
blocks with axial air passages, to which the friction brake 
block liners are riveted to insulate the rubber gland fron 
the heat generated at the friction surface (cf. Reference 22). 

Rubis was able to cbtain from the Falk Corporation the 

torque expressions for the clutches.(11) The torque capacity 


was found to be directly related to the number of brake liners 








PALK BARINE ALUYLia CLUTCH 


Fige (3-1) 
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3e CLUTCHES AND HYDRODYNAMIC TORQUE CONVERTERS (Cont'd) 


362 Falk Marine Airflex Clutch (Cont'd) 


and the net air pressure used to engage the clutch. For a 
35~inch-diameter clutch with eight liners, the maximum clutch 
torque capability, Qe» is given by the expression 

Q, = 2275 Pia, (lb-ft) 


p 


net = het air pressure activating the clutch (psi) 


One of the important features to consider in this clutch 

is that the inward expansion of the clutch glands during 
engagement is opposed by the centrifugal force of the rotating 
clutch assembly. The retarding force has been transformed 

to an equivalent air pressure, pc, and is expressed by the 


relation 


pe = 5 + (5.9 x 107?)g” (psi) 


Ng = speed of the frame housing the clutch glands (pm) 
Rubis states that typical inflation and deflation rates are 
5 psi per second and 30 psi per second respectively; and a 
typical supply pressure is 150 psi. 
SQ500 Synchroclutch (Dental Clutch Type) 

The SQ500 Synchroclutch is manufactured by the Philadelphia 
Gear Corporation. [Literature on the description and operation 
was furnished by the manufacturer and is included in Appendix C. 
This particular clutch is recommended for gas turbines in the 
25000 HP at 3600 rpm range (e.¢. FI4A-2 & IN2500). Fora 
typical application, the torque history for the synchronizing 


(friction) clutch may be approximated by: 





ad 


3. CLUTCHES AND HYDRODYNAMIC TORQUE COIMVERTERS (Cont'd) 


303 


SQ 500 Synchroclutch (Dantal Clutch Type) (Cont'd) 


_ t(qeo7t/26) 


torque at time t (ft~lb) 


ul 


F 


T 
2 
t time (sec) 


r 


Both r and the time constant may be changed to suit specific 


11,600 (ft-lb) 


if 


operating conditions.(23) Other information end details on 
this clutch are included in Appendix C. 

Another clutch which operates in a similar manner to the 
SQ500 is the BLH-Dynetic synchronizing clutch manufactured 
by Koppers Co, Inc. Although no torque-time history expression 
was provided by the Soriganuaaaie information was provided 
on the operation and design of the clutch. This is included 
in Appendix C for information purposes only. 
Hydraulic Transmissions 

One of the problems inherent in a plant driven by a diesel 
or gas turbine prime mover with a fixed pitch propeller that 
is not found in a comparable steam plant is that the former 
prime movers have minimum speeds at which they can be operated. 
The Coast Guard Ship Hamilton (WPG-715), for example, is powered 
Gy eee FRteA-6 gas turbines and two Fairbanks Morse Model 
387D-1/8 turbo-supercharged opposedepiston marine diesel engines 
driving two controllable pitch propellers. With the propeller 
at full pitch, and the diesel at idle, the vessel will attain 


a speed of approximately 12 knots, and 18 knots with full pitch 
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3.4 Hydraulic Transmissions 





and turbine at idle.(24) Richardson describes a method 

whereby the Falk clutch can be partially inflated thereby 
permitting continuous peCret onnee maneuvering at shaft speed 
below that corresponding to minimum turbine (or diesel) idle.(25) 
For some applications, however, the heat rate may bs prohibitive 
for continually doing this, and its effect on the life of the 
Clutch is not clear. 

In Reference (26), Rowen describes a hydraulic trans- 
mission system that permits propeller speeds in the low rpms 
with a diesel or gas turbine. The system has a hydraulic 
reversing gear which incorporates a fluid coupling for ahead 
rotation of the propeller and a reversing torave converter for 
astern operation. In addition, a positive drive clutch is 
included for normal oneration in the ahead direction at sea. 
Figure 3-2 is a cross section of the hydraulic transmission 
system. The torque characteristics on a “per unit" basis of 
the two hydraulic elements are presented in figure 3-3.* 

Speed ratio, as defined by Rowen, is the ratio of the output 
speed to the input speed (positive for ahead rotation, negative 
for astern rotation). He assumed that the hydraulic element 
efficiencies were constant over the speed ratio range, and 
defined the output torque by the relation, 

Ze 


Ty i. Qor = Pom 8 


* Per Unit" is defined as the actual. value divided by the design point value. 
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3e4 Hydraulic Transmissions (Cont'd) 
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iis = converter or coupling output torque 
= oil flow to converter or coupling 

Nom = input speed of prime Feuer 

F =a function of the speed ratio (Fig. 3-3) 


Figure 3-!1 is a block diagram of the hydraulic transmission 
simulation used by rowen Sto nunbers were given by him for 
the various parameters involved in the simulation. 
Discussion 

The Falk clutch and the Philadelphia Gear synchroclutch 
were selected for inclusion in the simulation because of their 
simplicity and because of their wide spread use in ship proe 
pulsion trains. Since no expression for clutch disengagement 
was provided by Philadelphia Gear, it was assumed that the 
clutch torque decreased as an exponential decay with a time 


constant of 5 seconds. 
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4.1 Introduction 

The propeller simulation proved to be the most challenging 
aspect of this study. An excellent study of propellers under 
conditions of reversing was meen by Miniovich and is listed 
as Reference (28). Miniovich questions the validity of using 
propeller characteristics measured under equilibrium conditions 
in the simulation of backing and reversing of ships. To answer 
this, he performed reversing experiments with models in a towing 
tank and measured the propeller characteristics. These char- 
acteristics were found to be in agreement with those measured 
in equilibrium or quasi-stationary tests. The only notable 
discrepancy appeared to lis in the zero thrust region. From 
his equations of motion, Miniovich also caiculated the values 
of basic variables that characterize the reversing process 
(e.ge, time reavired to stop the ship, time reauired to stop 
the propeller, and maximum values of thrust and torque) and 
found these to differ from experimental data by a maximum 
of 10 percent. It thus appears that one can use ths steady 
state propeller characteristics in a simulation and get a 
reasonable accuracy. 

4.2 Propeller Characteristics 

‘The traditional definition of the propeller advance 

coefficient is the familiar expression, 


x 2A 


ND 


V, = propeller speed of advance (ft/sec) 
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4. PROPELLER MODELS (Cont'd 


4.2 Propeller Characteristics (Cont'd) 


N = propeller rate of rotation (rev/sec) 


D propeller diameter (ft) 

This coefficient can be derived — a dimensional analysis 

for propellers (cf. Reference 29). It is obvious that relating 
fixed pitch propeller coefficients to J might prove to bs 
unwieldy for ship propulsion studies since J becomes unbounded 
as the shaft goes from positive rotation (providing thrust 

to move the ship ahead) to negative rotation (prceviding thrust 


to move the ship astern). Miniovich solved this dilemma by 


defining the propelier characteristics in terms of both J 
ND 
Vas The result is two graphs of 


propeller characteristics instead of ones when either one 


and its inverse, Vize, 


of the coefficients tends to large values he changes to the 
alternate definition. 

At the Ninth International Towing Tank Conference (1957), 
Stephanson demonstrated how the J and the 1/J propeller 
characteristics could be combined with a new term, called 


the modified advance coefficient defined by, 


a = ote 
Cv, Neo 


where all the terms are defined as before.(30) (Miniovich 
ascribes a similar expression to Kalmakov in 1940). It can 
be seen that J* will not become unbounded as N or Vy approach 


Zero valuese 
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“Wn, PROPELLER MODELS (Cont'd ) 
4.2 Propeller Characteristics (Cont'd) 


In Reference (31) Baker and Patterson discuss Miniovich's 
work and present his propeller characteristics using the modie 
fied advance coefficient. Other forms of advance coefficients 
are offered by the two authors to make the task of simulating 
a propeller easier than as performed by Miniovich. 


Modeling the Waveningen ReScrew Series 





In 1969, Van Lammeren, Van Manen, and Oosterveld presented 
a new techniaue of reproducing the propeller characteristics 
of feiaecewasen B-Screw Series.(32) Ths nondimensionalized 
thrust and torque ssskaueeniee’ presented in the forn of 
a Fourler series. The series a as a function of the 


hydrodynamic pltch angle, Bs defined by, 


V 
p = arctan ( “A 
0.7 Of ND 
where Vails and D are as previously defined. The nondimene 


sionalized thrust coefficient, Cs and torque coefficient, 


C,*s are then given by 


Ct = M 
T ae , 
Koo. CAM(K) cos BK + B(K) sin AK 
loc*= Mf 
=o A(K) cos BK + B(K) sin AK 


Once the coefficients are obtained, thrust, T, in pounds, 


and torque, Q, in ftelbs can be calculated using, 


Se 4 8p Ke + (0.71 D aes ef 


4. PROPELIER MODELS (Cont'd) 
4.3 Modeling the Wageningen B-Screw Series (Cont'd ) 


~ 7 a 2 | GAT 3 
Q = ctf pL V, rn (0.7 NT D Few 
where 
ie = fluid density 
V,; N; Das previously defined 


A 


The factor "10" in the expression for noncdimensionalized 
torque coefficient apparently through ommission, does not 
appear in two different published versions of Reference (32). 
Simulations were performed without this factor and resulted 
in excess propeller torques (10 times that predicted from 
hand calculations). Figure (40) of Reference (32) was 
reproduced by the IPM 1130 computer without the factor of 

10 in the expression. Eaten of Figure (40) revealed 


i 
being introduced to avoid overlapping of curves. It was thus 


that the ordinate is labeled C,* and -10 Cys the minus sign 


deduced that the expression for Sa should include a factor of 10. 
The values of the coefficients A(K) and B(K) for a series 

of 20 terms, i.e. K = 1-20, are presented in Reference (32). 

These coefficients are presented for a wide range of B-Screw 

series. Figure 4-1 illustrates how closely a series of 20 

terms matches the measured characteristics of a B 4-70 screw 

with a pitch ratio of 1.0. Figure 4-2, furthermore, illustrates 

how well a series of 10 terms approximates the measured chare- 


acteristics for the same propeller. 
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4, PROPELLER MODELS (Cont'd) 


4.3 


Lh 


Modeling the Wageningen B-Screw Series (Cont*d) 


It was probably noticed that the definition of 6 does 
not avert the problems previously associated with the tradi- 
tional definition of J, i1.c., dividing when n=o. It would 


seem more appropriate to use 


vA 
« 2 Cy 
arcsine (v,* + (Cad WY IID) yz 


I 


if 


p 


= arccosine (V,~ + (0.7 4 ND)~) 


4 
2 


instead of the arctangrent definition of ee The IBM 
scientific subroutine library on tne IBM 12130 digital comouter 
that was used did not have subroutines for calculating arcccosine 
or arcsine. The arctangent subroutine had to ce used. It 
should be also noted that the IPM subroutine for arctangent 
generates values only between Z WF to + [GC iS the 
Z Ze 

writer's opinion that if a numerical overflow should occur 
in the argument for tne arctangent subroutine, that this sub. 
routine will yield either ™/2 oy = 1/2 depending on the signa 
of the overflow. No numerical overflow problems were serenteneed 
using the arctangent subroutine 
Controllable Reversible Pitch Propellers (CRPP) 

A simple model of a CRPP was presented by Rowen in 
Reference (26). The "per unit" equation for thrust is given 
by 

Thrust = (pitch x rpn)* 





4%. PROPELLER MODSLS (Cont'd) 

4.4 Controllable Reversible Pitch Propellers (CRPP) (Cont'd) 
and the "per unit" equation for torque is given by 
Torque = (0.3 x (ram? | 7 (0.7 + (eae x (piteh)?] 
In Reference (33), Read presents a stndy of gas turbines 
driving a CRPP, and includes the torque and thrust characteristics 
for the propeller. Apart from the characteristics, all that 
is known about the propeller is that it has 4% blades and has 
@ polar moment of inertia of 7050 Lbeftesec”, Wendel and 
Dunne apparently found a small error in the values of medified 
torque coefficient piven by Read and applied a correction to 
the torque characteristic.(34) The modified advance coeffi- 
cient, J', is as previously defined, viz., 

J§ = Vy aes 


(ee + N“D~)? 


The modified thrust coefficient is defined by 
i ; 
a Fo, + ND’) 


and the modified torque coefficient is defined by 


K = al Q 
i r,* + N°D*) 
where T= thrust in lbs 
Q = torque ft-lbs 
fs Vy» N, D are as previously defined 
(Note the absence of 77 in the denominator in the expressions 


for Kn" and K's ff 4s included in the similar expressions 
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4. PROPELLER MODELS 





4.4 Controllable Reversible Pitch Propellers 





given by Van Lammeren. Baker and Patterson discuss the 
w question" in Reference (31)). 

Figures (4-3) and (4-4) are the CRPP characteristics, and 
include the Wendel-Dunne correction. 

4.5 Pitch Control System Model 

The only analytical description of a CRPP pitch control 
system found was that described by Read in Reference (33). 
Pitch, Ry and angular rotation,§, are related by the expression 


R = Tan. @ 
00455 


The pitch-changing cylinder has a total capacity of 40 gallons, 
and the blade angle varies very nearly linearly with ram stroke, 


such that 


dé = % FE (Rad/sec) 
‘ole 57 3 18 


Fo = oil flow in gal/sec 
By manipulating the two equations, one can arrive at the 
expression 


dp = + 0.0528 + 0.0109 RF 
at 


There are two pumps supplying the oil flow, F One is 


O° 
motor driven and the other is shaft driven. The motor 

driven pump has a maximum capacity of 60 gal/min and the 
shaft driven pump has a maximum capacity of 100 gal/min. 


The pumps are pressure controlled such that any alteration 


in pitch demand causes the motor-driven pump to provide 60 
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4, PROPELLER MODELS (Cont'd) 


4.5 Pitch Control System Model (Cont'd) 


4.6 


gal/min and the shaft-driven pump to provide its maximum 
output at the shaft speed in question. The result is that 


during a pitch alteration, 


Fy = (1 + 0.435 m1) (gal/sec) 
N = shaft speed (rev/sec) 
Discussion 


fhe fixed pitch propeller model selected was the Fourier 
representation of propeller charactcristics as presented in 
Reference (32). Only the first ten coefficients were used 
in the representation. The resistance characteristics of the 
vessel, and the propeller diameter are two important factors | 
in selecting the proper propeller. Since the characteristics 
for the B-4 series with an expanded area ratio of .7 were 
readily available, it was decidei to use this type of propeller. 
A plot of K n/a” indicated a pitch-diameter ratio of 1.0 would 
give the maximum efficiency; hence this propeller was chosen. 
Read's controllable reversible pitch propeller (CRPP) 
was chosen. The propeller characteristics as modified in 
Reference (34%) were redrawn for constant values of modified 
advanced coefficient (J*). The independent variable was 
taken as pitch diameter (PD) ratio and the dependent variables 
were thrust and torque coefficients. Figure (4-3) shows the - 
thrust coefficient vs. PD ratio for J' = -0.2. This curve 
is fairly typical of the others. The linear approximations 


are indicated. Figure (4-4) shows the torque coefficients 
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4, PROPELLER MODELS (Cont'd) 
4.6 Discussion (Cont'd) 
vs PD ratio for J' = -0.2 with the linear approximations 
used. This curve is also fairly typical of the other 
curves of torque coefficients for constant J' and varying 
PD ratio. 
Read's pumping system as described earlier was also 


used in conjunction with the CRPP. 





dy 


5 SHIP CHARACTERISTICS 
5el Wake Fraction and Thrust Deduction Factor 
Wake fraction, w, determines the velocity the propeller 


139 advancing at in relation to the ship velocity, 


Vi, = (lew) V 

Vy = propeller speed of advance 
V = ship speed 

w = wake fraction 


Thrust deduction factor, t, indicates how much propeller 


thrust is reduced due to hull interaction, 


Sec. (1-t) | 

ow = Open water propeller thrust 
t = actual thrust 

t = thrust deduction factor 


(For a physical explanation of "t" and wt, and other useful 
relationships involving these ae terns, see Reference (29), 
Chapter VII, Section 13). 

The implication in the two above definitions is that 
they are steady state quantities. Thus at equilibrium, the 
actual thrust equals ship resistance. Since t<<1l., the 
definition shows how much more thrust is required over resis- 
tance to keep ship speed constant. 

Wake fraction and thrust deduction factor are not, as 
one would expect, constant over a range of speed. The variation 
in t and w depend on the degree of interaction between propeller 


and hull. Harvald performed several towing tank tests to 





5 


5e SHIP CHARACTERISTICS (Cont'd) 
51 Wake Fraction and Thrust Deduction Factor (Cont'd) 

measure "t" and “wt under transient conditions for two 
different hull shapes.(35) He performed his towing tank 
tests under conditions of constant model ship velocity and 
varying shaft rpm and vice versa. Although his work pro- 
vided good insight into the behavior of "t" and 'w't under these 
two conditions, it was difficult to assess how his results 
coula be applied for conditions of varying shaft rpm and Pe. 
velocity. How did others treat this phenomena? In the solue 
tion of his backing equations, Miniovich assumed both factors 
to be equal to zero.(36) RBubis assumed varying values of 
ute and wu" from two curves of bell shaped form.(37) Read 
considered wake fraction to be negligible, but formulated 
linear relations for thrust deduction factor for positive 
and negative values of propeller thrust.(38). Lewis, Lecount, 
and Scoville assumed thrust deduction factor and wake fraction 
to be a function of vessel speed only; for negative speeds 
thrust deductions and wake fractions were neglected.(39) 
Thus, it can be seen that there is no general agreement on 
the treatment of these two factors. The feeling appears to 
be that these two factors are so small that any errors in 
their treatment will have negligible effects on the total 
simulation. For simulations performed here, constant alae 
of wake fraction (w=0.04) and thrust deduction factor (t=0.02) 


were assumed. 





5. SHIP CHARACTERISTICS (Cont'd) 


5«2 Added Mass of Ship and Propeller 
As with thrust deduction and wake fraction, the treat- 


ment of added mass of ship and propeller is also varied. 

The added mass of the ship is usually taken to be between 
82-10% of the displacement. This study used 10%. Lewis, 
Lecount, and Scoville used a different approach to account 
for added mass as will be discussed later. The added mass 
of a fixed pitch propeller is usually assumed to be about 
25% of the inertia. (cf.,@efe, Reference (40)) Twenty-five 
percent was assumed here. There is nothing fixed about this 
figure as Burrill and Rebinson found the percentage additicn 
should be as high as 46% for one particular propeller.(41) 


For his CRPP, Read proposed in Reference (33) the expression, 


E = 325 +2050 (Ry 
I = added mass due to entrained water (1b-ft-sec*) 
R = pitch diameter ratio 


This expression was used in the simulation with the B&W 
7K98FF low speed diesel ~ CRPP combination. 
503 Ship Resistance 

Accurate representation of ship resistance is, of course, 
vital to the proper selection of propeller(s). One of the 
basic assumptions of the popular "Propeller Law* (cf. References 
(14) & (15) ) is that resistance is proportional to the square 
of velocity for ships where the wave drag is small. For those 
cases where it is desired to simulate transients on a parti- 


Cular ship, it is, of course, more accurate to use that 
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5. SHIP CHARACTERISTICS (Cont'd) 


Ship Resistance (Cont'd 

vessel's characteristics as determined from towing tank tests. 
Lewis, Lecount, and Scoville included acceleration effects 

in their expression for ship resistance, 


a dv 
Ry Flees sie — 


M mass of ship 


Ne = ship speed 
ee ae (V.) = steady-state resistance 
Can = “Nadded mass" coefficient 

No values for C,. were provided by the authors. 


AD 
For simulations performed here, resistance was assumed 


to be of the forn, 


R =10%.8V* (lbs) 


ship velocity (ft/sec) 


tt 


V 
Sheft Friction 
Read used the expression 


2890 nt*?3 (Leb) 


Qe 
N 


propeller speed (assumed to be rev/sec) 
to represent shaft friction.(33) Rubis, on the other hand, 
represented shaft friction characteristics in tne form of a 
"table-look=up" in his computer simulations. These characteris- 
co iae approximated in simulations performed here by the 
expressions 

Q- = 2500 (ft-lb) 0 </Nn/ < 25 

Qp = 100/n/ (ft-lb) /n/ > 25 


N = propeller shaft rpm 
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6.1 Introduction 


642 


6.3 


The following maneuvers were performed with the drive 
trains as indicated: | 

FIWA-2 Gas Turbine - Fixed pitch propeller 

le Coastdown from 30 knots 

2e Accelerating from 5 knots 


3. Crashback from 30 knots 


L7Z2500-A Ga Gas Turbine « controllable reversible pitch 
propeller — 


le. Accelerating from 15 knots 
Le Crashback from 30 knots 


Mirrless KV Major 12 «= Fixed Pitch propeller 





le. Crashback from 9 knots 


Burmeister & Wain 7K9SFF - Controllable reversible pitch 
propeller 


1. Crashback from 12 knots 
Coastdown from 30 knots (FT4A=2 

At zero time, fuel rate was decreased from 9500 lb/hr to 
1660 lb/hr exponentially with a lesecond time constant. 
Results of this simulation are presented in figures (6-1) - 
(61+). 
Accelerating from 5 knots (FT4A-]=2 

At zero time, fuel rate was increased exponentially with 
a lesecond time constant from 2500 LB/hr to 12000 LB/hr. 
Figures (6-5) - (6-11) are plots of the results. A copy of 


the computer program used is included in Appendix F. 
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6. RESULTS (Cont'd 
6. Crashback from 30 knots (FTA-2) 
This maneuver was accomplished in three phases: 
1. Coastdown* - opening clutch 
Ze engaging reverse clutch 
3- stopping the ship 
Results of this simulation are plotted in figures (6-12) ~< 
(6-17). Time for forward clutch disengagement was 4.8 seconds 
The reverse clutch was engaged at 11.2 seconds and the ship 
stopping time was 43.4 seconds. Head reach was 1229 ft. 
A copy of the computer program used in phase 3 is included 
in Appendix F. 
6.5 TM2500~/ ) 


Accelerating from 15 knots 





A description of the control. strategy and the controller 
used is in Appendix G Initially the ship is at 15 lmots 
with a pitch ratio of 1.0 on the CRPP. Pitch was increased 
to 1.6 in 8 seconds. At that instant the controller acted 
to increase shaft rpm to 216. Resvlts of this simulation 
are shown in figures (6-18) = (6-20). 

6.6 Crashback from 30 knots (LM2500-A) 

The problem to consider here was whether or not the power 
turbine would overspeed while the CRPP goes through a zero 
value of torque coefficient. This ship is initially at a 


speed of about 30 knots with full pitch ratio (1.6). The 


*Coastdown as used hereafter implies setting fuel rate to idle. 





6. RESULTS (Cont'd) 


6 6 


Gu? 


Crashback from 30 knots (LM2500~A Cont'd 


first simulation that was tried used the following sequence 
of events: 

ae reduce fuel to idle 

b. decrease pitch to -1.2 

c. increase fuel to 9900 lb/hr 
In performing this simulation, it was observed that the pitch - 
changing operation took a rather lengthy time. The end result 
would have been a time to stop in excess of 60 seconds. The 
simulation was stopped and a new sequence of events was 
attempted, 


&e Simultaneously decrease fuel to idle and decresse 
pitch 


be when pitch equals -.8, increase fuel to 9900 
Using this sequence of events, 'a stopping time of 44.3 
seconds resulted. Results of this simulation are plotted 
in figures (6-24: - (6624). A copy of the computer program 
used in sequence (a) is included in Appendix G. 

Crashback from 9 knots (KV Major 12) 

The same sequence of events described in section 6.4 
were used here, ViZe, 

1. coastdown = opening clutch 

2e engaging reverse clutch 

3. stopping the ship 
Forward clutch disengagement occurred at ll seconds, and 


6.2 seconds later, the reverse clutch engaged. The ship 
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6. RESULTS (Cont'd) 
One Crashback from 9 knots (KV Major 12) (Cont'd) 


6.8 


was stopped in 95.2 seconds with a head reach of 1208.7 feet. 
Results of this computer simulation are plotted in figures 
(6-25) ~ (6-27). Copies of the programs used in all three 
sequences are included in Appendix H. 
Crashback from 12 knots (B & W 7K98FF) 
This simulation successfully illustrated the increase 
in prime mover speed as the CRPP went through a zero torque 
coefficient region. Tne sequence of cvents employed here were: 
1. simultaneously decrease fuel to idle and pitch ratio 


Ze when pitch ratio equaled «-.2, increasa fuel to stop 
the ship 


Tne initial simulation started from a pitch ratio of 1.0 and 
decreased pitch to -1.2. At this pitch ratio, the propeller 
torques were so high that the equivalent. of diesel stalling 
occurred. The change in pitch ratio was modified to go fron 
+1. to -l. This proved to be satisfactory. The time of pitch 
change was 34.5 seconds. Ship stopping time was 90.5 seconds 
with a head reach 629.8 ft. The results of this simulation 
are plotted in figures (6-28) = (6-31). A copy of the pro- 


gram used in sequence (2) is included in Appendix I. 
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CONCLUSIONS 


Simple models were found to exist that enable one to model 
ship propulsion drive trains. Some of the models had to be refined 
by the addition of time lags, delays, etc. ‘he combination of 
drive train components must be compatible with the characteristics 
of each drive train member. The simulations illustrated the 
importance of proper sequencing of propulsion plant control para- 
meterse Conditions of overspeeding, excess torques, and high 
energy Gissipation rates in clutches can be easily simulated. 

These simulations can prove to be immensely useful to the control 
engineer in Hel pineaiae synthesize a controller that will keep a 


system within its region of safe operation. 
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A study should be taken to determine what is the optimal 
fuel flow Getting) - pitch ratio combination for a Pe eiar 
prime mover and the CRPP presented here. 

A model of reversing the B & W 7K93FF should be developed 
and a crashback simulation of the 7K98FF with a fixed pitch 
propeller should be performed. 

A model controller to keep shaft rpm nearly constant while 
changing pitch ratio on the CRPP could be synthesized for 
each one of the prime movers. 

A torque converter model. should be developed to be used with 
the prime movers developed here, and maneuvers should be 


simulated. 
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APPENDIX A 


Gas Turbine Models 


Lurye Model 
The following model of a marine gas turbine was proposed by 
Larye in Reference (10) in the form of linearized differential 


equations: 


Turbocompreéssor 


ae _ Se 
(Ts +1) 4% = Ki AG * KA 


Free Turvine 


=7, 
(T, Ss + 1) Aw, = K,,4 G2 + K,,AMw =~ by A fh 














qT = vA eoee Unrbocompressor 
_ IM ») e/4,. time constant 
J iy, Ite, 
A/1, 

44 = o} coc ampauification coefficienc 
ol ty a, of turboccnipressor by fuel 
dt, Oe, 

SFb 
Kyo = g kis __. eee load coefficient by the rate 
fs oe J, of rotation of free turbine 
a 
o) Ce, ow, 

poe = A, ooee time constant of free 

é _ AM, P oF turbine 
Ju, Dp ce, 

ees 

Kaq = od G, occ AUPLilication coefiicieny 

eh d/4, M7, of free turbine by fuel 
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(Cont'd) 
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APP&NDIX A (Cont'd) 











Mz 
op eeee Load coefficient by propeller 
is aM, y Le, pitch 
ot, Aas, 
A, 
jan esos amplification coefficient by 
rate of compressor rotation 

SAM, _ 4, mp 
Je, oa 


turbocompressor rate of rotation 
free turbine rate of rotation 
fuel consumption 

propeller pitch 


time lag in combustion chamber 


“ty (wy, ,G) torque moment of turbocompressor 


ee | 
oY) 


w ' ; 

1, Ws, Fo) torgue moment of free turbine 
turbocompressor rotational inertia 
free turbine rotational inertia 


f.. (w, , 9 ) resistance moment of propeller 


f (4, WwW.) resistance moment of compressor 
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APPENDIX A (cont'd) 


An Extension of the Saravanzinnttoo Model 
Using the technique proposed by Saravanamuttoo in References 
(7) & (8), the model can be extended to the case where 2-spool gas 


generatorefree turbine combination: 


q, 


HP Rotor equation CZ, 
a he Be a 
Az, yy Co pr (A= ip) - oes t HPC pe GE) ae 
LP Rotor equation Mee 
%r, Mp7 Ger Cie ~ Ze) — <p “W476 CF pe CE-T) = = aes 





ha Tucbine eauation 
Croc. Mpp Cp pp (Mg > Ty) 
LP-HP Intercompressor Volune 


RL, 6) 
we \J — 2 
Wupc 7 "ppp = —& —-& 


Vo abt 
Combustor volume 
: . ee ay 
Wupr 7 “upc 7 @y oo Eh 
Ky Yta 


Interturbine Volume (LP-HP) 
. wees 
roe ee 
Ua ofA 

Interturbine Volume (HP-FT) 


ee A 
Wuep > “pr > Se He 
7 at 


W 
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| APPENDIX A (Cont'd) 


Le (Cont'd) 


Combustor 
eae big AVA Cp, Wyse Ap 
op, (My + Uap) 
Nomenclature 
W = mass flow rate 
T = temperature 
c= specific heat at constant pressure 


#,= turbine over-all efficiency 

%, = compressor over-all efficiency 
K = fas constant 

HV= fuel heating value 

V¥ = volume 

NW = rotational speed 

P = pressure 

t = time 

I = polar moment of inertia 


Station Numbering 


1 = LP Compressor inlet 
2 = HP Compressor inlet, 
3 = combustor inlet 


4 = HP turbine inlet (combustor cutlet) 
5 = LP turbine inlet 
6 = FT inlet 


7 = FT ontlet 











Ze 


(Cont'd) 


subscripts 


LPC 
HPC 
FB 
HPT 
LPT 
FT 
HP 


i 


tt 
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APPENDIX A (Cont'd) 


lew pressure compressor 
high pressure coupressor 
Combustor 

high pressure turbine 
low pressure turbine 
free turbine 

HPCeHFT rotor 


LFPC-LPT roter 
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APPENDIX A (Cont'd) 


FTMA}-2 & LM2500-A Marine Gas Turbines 





Figure (A-l) is the pover«speed curve for constant SFC valucs. 
The torque-speed characteristics shovm in figure (A-2) can be 
derived from this curve. Figure (A-3) and (A-4) are the corres~ 
ponding curves for the 1M2500-4A. Figures (A-2) and (A-4) show 
examples of the linear approximations for constant fuel rate 
values used to approximate the gas turbines‘ characteristics. For 
intermediate values of fuel rate, linear interpolation was uscd. 
The subroutines used for the FRA-2 and the LM2500-A are included 
in the rear of this appendix. 

A sketch of the PRM AW2 with a JFC 25828 fuel controller is 
shown in Figure (A-5). <A controller similar to the JFC 25-28 is 
the JFC 25-7 shown in figure (A-6). Tais controller is typical 


of the non-linear systems engineers are often asked to model. 
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l 057-637 


SUBROUTINE LM250 (RPM,WFE, TORQ, A08,A]69A21A3:A4,A5, AG AT,A99) 
THIS SUBROUTINE CALCULATES FREE TURBINE TORQUE AS A FUNCTION OF 
RPM AND FUEL FLOW 

MINIMUM FUEL RATE IS 800 LP/HR 

MAXIMUM FUEL RATE ALLOWABLE IS 9900 LB/HR 
ENTERING ARGUMENTS FOR THIS SUBROUTINE ARE FUEL RATE (WF) 
AND PQWER TURBINE RPM (RPM) OUTPUT IS POWER TURBINE TORQUE IN 
FT-LBS 

SEARCH TO BRACKET WF AND THEN PROCEED TO FUFL SUBROUTINES 
Mepdr «EO. 800.) GOTOLS5 

Mmetnwr .L&. 1600.) GOTO1L6G 

Mere .LE. 2000.) GOTOL? 

IF(WF .LE. 3000.)GOTO18 

Meiewr .LE~ 4000.)GCIOLO 

Meer oLE. 5000.) GOTO 11 

TFCWF .LE.6000.)GOTOL2 

Meenr LE. 7000.)GOTO1L3 

IF(WF .LE. 9900.)GOTO14 
TORQ=-1.965#RPM+5500. © 

RETURN 

CALL LWFOS8 (2 PI’, A08) 

CALL L¥YF16 (RPM,A1L6) 
TORQ=A08+(WF-£00.)*.00125%(A16-A08) 
RETURN 

CALL LWF2 (RPM,A2) 

CALE LWF1G (RPM,A16) 
TORQ=ALG+(WF-1600.)%.0025%(A2-Al16) 

RETURN 

CALL LWF3 (RPM,A3) 

eel LWF2 (RPM,AZ2) 

TORQ=A2+i(WF-2000.)% -001%(A3-A2) 

RETURN 

CALL LWF4 (RPM,A4) 

CALL LWF3 (RPM,A3) 
TORQ=A34+(WF-3000.)#. 0018 (A4-A3) 

RETURN 

CALL LWFEF5 (RPM,A5) 

CALL LWF4 (RPM,AG) 
TORQ=A4+(WE-4000.)#-00L#(A5—-AG) 

RETURN 

CALL LWF6 (RPM,AG6) 

CALL LWFS (RPM,AS) 
TORQ=A54+(WF-5000.)#.001#(A6—-A5) 

RETURN 

CALL LWFY (2PM,AT7) 

CALL LWF& (RPM,AG) 
TORQ=A6+(WF-6000.)#.001%(AT—-AG) 

RETURN 

CALL LWF99 (RPM,AQ9) 

CALL LWF7 (RPM,AT) 

TORQ=AT7+(WF- 7000.) #*.0003448#(A99- AT) 
RETURN 

ENO 

SUBROUTINE LWFO8 (RPM,AOB) 
AO08=-1.965#RPM4+5500. 

RETURN 





PAGE 
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ENO 

SUBROUTINE LWF16 (RPM,A16) 
A1l6=-2.94#RPM+12500. 
RETURN 

END 

SUBROUTINE LWF2 (RPM,A2]) 
A2=-3.48#2PM+16150. 

mE TURN 

END 

SUBROUTINE LWFEF3 (RPM,A3) 
A3=-4.77#RPM4+23250. 
RETURN 

EWD 

SUBROUTINE LWF4 (RPM,A4) 
A4=-5.933#RPM+30800. 
RelUAN 

END 
SUBROUTINE LWEFS (RPM,ASD) 
A5=-7.133#RPM+329500. 
RETURN 

END 

SeeeRUUTINE LWEF6 (RPM,AG6) 
No=-/.5#RPM+47500. 
RETURN 

END 

SUBROUTINE LWFE7 (RPM,A7) 
A7=—-8.#RPM4+55300. 

RETURN 

END 

SUBROUTINE LWF99 (RPM,A99) 
A99=-10.15#2PM+73600. 
RETURN 

END 
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APPENDIX B 


Diesel Particulars 


Burmeister and Wain 7K98FF 

Test-hed results of the B &W 7K98FF are shown in Figure 
(B-1). From this figure, it was possible to derive the ideal 
diesel power-speed curves for constant fuel settings; this is 
shown in figure (B-2). The torque-speed characteristics, figure 
(B-3), are in turn derived from figure (B-2). The "FF" at the 
end of the designation of this enzine denotes that this diesel is 
intended for direct drive, i.e., no reduction gears. Other 
technical data is presented below: 

continuous service rating 24500 metric BHP @ 100 rpm 

maximum continuous service rating 26600 metric BHP © 103 rpm 

Length 45* 2" | 

Breadth 16' 2" 

Height 43° 2 

The computer program used to model the diesel is included 
at the end of this appendix. 
Mirrless KV Major 12 

A cutaway drawing and technical specifications of this engine 
is shown in ek (B-4). The "12" at the end of the designation 
indicates the number of cylinders. Figure (B-5) shows the fuel 
map for this engine. Using this fuel map and the elementary 
diesel equation discussed in section 2, it was possible to obtain 
the approximately Linear torque-fuel setting characteristics shown 


in figure (B-6). The relationship obtained from figure (B-6) is 
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APPENDIX B (Cont'd) 





Ze (Cont'd) 
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APPENDIX C 


Clutches 


Figures (C-1) and (C-2), which were provided by the Philadelphia 
Gear Company, describe the operation of quill shaft version of the firm's 
synchcroclutch. Additional useful information provided by Philadelphia 


Gear Is; 


2 


MHertia ac se See sees cee oa —1 4 lb-ft, 


Output -110 Ib-ft 
torsional stiffness*eesesecces Approx. 270 million inwlb/Rad 
Figures (C-3) and (C4), which were provided by the Koppgers Co., 


describe the operation of the Blile-Dynetic clutches. 


* For information purposes only 
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ROTARY UNION 
OENTAL CLUTCH ASSEMBLY 


ENGAGING CYLINDER FOR DENTAL 
CLUTCH 
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. SHIFT COLLAR 
. PINS 

. SLEEVE 


INTERMEDIATE CEAR 


. DRIVE SPLINE 

. OUTPUT FLANGE 

. SELF-LOCK PISTON ASSY. 

. DISENCAGE CYLINDER-DENTAL 


CLUTCH 


. SPACE FOR OPTION. QUILL SHAFI 





BRAKE 
. FRICTION APPLY CYLINDER 


HYDRAULIC BALANCE & FRICTION 
RETURN PISTON 


Q. JOURNAL & THRUST BEARING 
P ASSEMBLY 


FRICTION CLUTCH 
COOLING OIL PORTS 9) 


Ww 


FIG. 2 CROSS SECTION OF TYPICAL SYNCHROCLUTCH (QUILL SHAFT VERSION) 


, me gy § ‘*: 

PDESCRIPTIC 
The new Philadelphia SYNCHROCLUTCH consists of a 
friction clutch in parallel with a positive drive dental clutch. 
The friction clutch synchronizes the driving and driven 
shafts, after which the dental cluich is engaged to provide 
a positive drive. The dental clutch features fully crowned, 
nitrided gears for maximum capacity for a given diameter. 
The friction clutch is of the multi-disc type with positive 
ol! circulation for long life and heat removal. The friction 
clutch has ample capacity to absorb thermal loads imposed 
by the synchronizing cycle—even very high loads encoun- 
tered when reversing large marine gears in the 40,000 
plus HP range. In addition, dead shaft pickup is possible 
in large. marine propulsion systems. This permits the prime 
mover to be run at speed in port for auxiliaries; such as 
generators, while retaining the ability to get underway 
instantly by clutching in the propeller. 


The dental clutch features a unique helical gear element 
which automatically aligns the teeth of the dental! clutch 
prior to engagement. This feature is illustrated in Fig. 3. 
The clutch nay be actuated hydraulically or pneumatically, 
using low pressure air (109 psi) or oi! (150 psi). The 
synchronizing action is controlled by a transistorized dif- 
ferential speed control using non-contacting magnotic 
pick-ups for speed indication and position readout. Once 
engaged, the dental clutch does not require power to 
remain in mesh, thus eliminating dependence on electrical 
or hydraulic power for continued torque transmission. 


MT AME @iels 


2 2 


* & 
b= fae £O 38 @ Paid 


External jacking is provided for emergency manua! engage- 
ment and disengagement of the dental clutch. 


Referring to Fig. 2, the friction clutch (a) is first engaged 
to force the driving and driven shafts into synchronism. 
The clutch is clamped by admitting oil (or air) to the piston 
through the rotary union (b). After the driving and driven 
shafts have reached a speed differential of approximately 
4%, the transistorized speed control signals for dente! 
clutch (c), engagement by admitting oil (or air) pressure 
to the engaging side of the dental clutch piston (d). The 
piston transmits the thrust through a heavy duty thrust 
bearing (¢c), to the rotating shift collar (f), which ts con- 
nected through large pins (g), to the sleeve (h), of the 
dentai clutch. The sleeve carries with it the intermediate 
gear (i), and meshes with the drive spline (j), on the output 
flange (k), in the manner iJlustrated in Fig. 2. A self-locking 
piston (I) on the dental clutch prevents piston back-off 
under torque reversal. This action is shown in Fig. 3. To 
disengage, pressure is admitted to the reverse side (m) of 
the dental clutch piston, which pulls the dental elements 
out of engagement as shown in Fig. 3. The friction clutch 
may then be de-energized to permit the driving and driven 
shafts to rotate relative to each other. Where desired, a 
quill shaft brake (n) may be incorporated into the clutch 
to stop the quill shaft from rotating, due to residual drag 
torque within the friction clutch. This feature requires no 
additiona! length. 


(C-1) 
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DENTAL CLUTCH TEETH. 
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ACCESSIBILITY. 


DIRECTION OF ROTATION. 


PERMISSIGLE AXIAL MOTION OF OUTPUT 
FLANGE CONNECTION. 


PERMISSIBLE RADIAL RUNOUT AT OUTPUT 
FLANGE. 

ELECTRICAL POWER REQUIREMENTS FOR 
CONTROLS. 

14. | CLUTCH ENGAGEABLE AT STANDSTILL. 


15. | EMERGENCY MANUAL CLUTCH OPERATION. 
OPERATION 
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4. SLIDING SLEEVE PUSHES INTERMEDIATE 
GEAR & DRIVE SPLINE TEETH INTO FULL 
MESH. 
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UP TO 100% OF RATED SPEED. 
OIL @ 150 PSI. (SPECIAL-AIR @ 100 PSI.) 

PGC. FSO] GRADE RCI. 
HARDENED STEEL. HIGH CARBON. GROUND. 


FORGED STEEL, HEAT TREATED AND 
DYNAMICALLY BALANCED. 


NITRIDED STEEL. _ 
CROWNED, CHAMFERED ON ENTRY SIDE. 
EXTERNAL—UNIT CONSTRUCTION, 


ALL CLUTCH MEMBERS REMOVABLE WITH- 
OUT REMOVING HUB OR QUILL SHAFT. 


Cw OR CCW. 
+1/8" STANDARD. 


.015 TIR ON 500 SIZE AND LARGER. 
010 TIR ONMiZ0 SIZE. 
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FUNDAMENTAL DESIGN 


All BLH-Dynetic Clutches are of the same basic 
design, modified as necessary to meet a wide vari- 
ety of application conditions. 

The forces developed during engagement are 
completely self-contained, and produce no thrust 
loads on the shafts or clutch housing. After engage- 
ment, power is transmitted through hardened and 
nitrided coupling teeth in positive mesh, with no 
loss due to slippage. Hydraulic pressure is required 
only during the actual engaging or disengaging 
cycle. The parts are held in their final positions by 
spring loaded detents, or mechanical locks, as 
necessary. 


DESCRIPTION 


OP bil 


MAJOR COMPONENTS 


Fa Operating Cylinder—to provide the axial force 
necesSary to engage or disengage the clutch. Actu- 
ation is usually hydraulic using lubricating oil as the 
fluid, but pneumatic operation is also avaiiable. 


By Friction Discs—to develop friction torque neces- 


sary to synchronize the shaft speeds. These discs 
are cooled by a continuous flow of lubricating oil. 


pel Transition Torque Control—spring loaded pins 
key with lugs in the load end hub to transmit torque 
during the transition from friction to gear tooth 
drive, and to align gear teeth prior to engagement. 


fy Gear Coupling—to transmit power after engage- 
ment is complete and accommodate shaft misalign- 
ment, if present. 


UY CLUTCi 


Fire (C-3) 
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Figure A 


POH 


Figure & 


Figure C 


Figure D 


OF BLU OD 


The heart of the synchronizing mechanism ts the 
cluster of synchronizing pins (SP) which cooperate 
with the holes in the power end hub (PEH) through 
which they project to control engagement as indi- 
cated by the diagrams at the left. 


As engagement begins the actuatirg piston moves 
the power end hub and its connected parts toward 
the load end hub. Through detent action, an initial 
friction torque is developed which swings the smail 
diameters of the synchronizing pins against thie 
sides of their respective holes in the power end hub, 
as in Figure A. 

Continued pressure brings the cam surfaces of 
the synchronizing pins into contact with those of 
the power end hub, thereby transmitting full piston 
pressure to the friction disc pack as in Figure B. 
The friction torque so developed acts to synchronize 
shaft speeds. The geometry of the assembly is such 
that the parts remain in the position shown as long 
as slipping occurs in the friction discs. 

Slipping stops at the friction discs when the shaft 
speeds synchronize. This unbalances the forces and 
permits the power end hub to advance over the syn- 
chronizing pins as in Figure C. At the same time, 
spring loaded pins are moved into spaces between 
lugs in the end of the load end hub, acting as loose 
keys to control torque until the coupling teeth 
engage. 

As the synchronizing pins become aligned with 
their mating holes as in Figure D, the cam surfaces 
separate, removing the axiai force on the friction 
disc pack. This permits the spring loaded pins to 
contact the sides of their lugs, maintaining syn- 
chronism of shaft speeds and aligning the coupling 
teeth for engagement as the piston completes its 
travel. 

With completion of engagement the hydraulic 
pressure is cut off, and the assembly transmits 
power as a simple gear-type coupling. 

During disengagement the piston moves the parts 
away from the load end hub, disengaging the spring- 
loaded pins from their lugs before the coupling 
teeth disengage, and securing the assembly in the 
final position to prevent dynamic unbalance. 


‘NGPIC CLUTCH 


Fige (CH4) 
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APPENDIX D 


Propeller Subroutines 


This appendix contains the computer programs used te model the 
p4-70 (P/D = le) propeller and Read's CRPP, both of which are discussed 
in Section 4. Included in this appendix is the program described in 


Section 5.4 that computes shaft friction. 
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BACE 1 | O37 027 


SUBROUTINE CRPP (J»CQ,CT,AP1,AP8,BP1,B8P8,AM8,ANL,BMES,BM1,AM4, BM4,A 
2142,BHN2;AM0,BHO,AP4,8P4,PD) 

CRPP COMPUTES MODIFIED THRUST) AND TORQUE COE Cle NS ce a 
VALUF OF MODIFIEDSADVANCE COEREITCIENT AND? It Cae eh aie oe 


B<OPELLOGR DATA 18 TAKEN GRO DYNAMIC. ANALY SESS a DiS tiie wie 
GF SHIP AND PROPULSION PLANT MANUVERING PERFORMANCE BY WENDEL 
moO DUNNE APPEARING IN PROCEEDINGS OF THE SECGND SHIP CONTE GE 
mec bEMS SYMPOSIUM. 


mMeetFITED ADVANCE COEFF.a...J=VA/SORT(VA®#24+(N#D) #42} 
MODIFIED THRUST COEFF...CT=T/(RHO#(De82) (VAR R2+(N#0) ##2) ) 
POOUIFIEFD TORQUE COEFF...~CQG=Q/(RHO#(D##3) e(VAH#24+(N*#D) # #2) ) 


RHO=DENSTTY OF WATER 

Pee ROUPELLOR DIAMETER 

feeee ROPELLOR SPEED OF ADVANCE 
fee RUPELLOR ROTATION RATE 
P=tHRUST 

Q=TORQUE 


Memes DIAMTER RATIOS ARE LIMITED BETWEEN +1.6 TO -1.2 
BerERING ARGUMENTS ARE MODIFIED ADVANCE CCEFFICIENT (J) AND 
geepch DIAMTEER RATIO (POD) 


ee I I Te Tk eC ee a ey) CP CY OY 


Be urTS ARE CQ AND CT 
REAL J 
Miao .lE.~.8)G0I1Q10 
ee LE .-24)G0TO11 
Mere lLE.--2)G0TO12 
mode LEO. )GOTO1L3 
IF(J.LE..4)G079N14 
Meds lLE.» ~-8)GOTOLS 
Meee PLUSIL (PD,AP1,BP1) 
Meee PLUSS (PD,;AP8,BP8) 
=n J-.6)4*5. 
MeeAPB&trs (APL-APB&) 
mr=BP8+F2(BP1-BPB) 
RETURN 

Memeercl MIONE (PD,AM!I,8M1) 
CALL MIN&S (PD,AM8,BM8) 
eo t) . KD. 
OI=AML +P ( ABS AMT) 
CT=BM)+F#(BM8-BM1 ) 
RETURN 

Me cALL MING (PD,AH4,BM4) 
CALL MIN8 (PD,AM8B;BM8B) 
= leas a) oe 
CQ=AMBtE# (AM4—Ai‘8 ) 
CT=BMB+F# (BM4—-BMB ) 
RUN 

cee ee Peay orc) 
CALL MINS (PD,AN4,BM4) 
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F=(J+.4)#5. 
CQ=AM4+4+F#(AM2-AM4 ) 
CT=BM44+F#(BM2-BMG ) 
RETURN 

BALL ZERO (PD,AMO,RMO) 
CALL MIN2 (PD,yAM2,BM2) 
Pedit.e2 ) #5. 

CQ=AM2+F# (AMO-AM2 ) 
Ol=BM2+F «+ ( BMO-BM2) 
RETURN 

CALL PLUS4 (PD,AP4,BP4) 
CALL ZERO (PD,AM0O,BMQ) 
feu t 7225 

CO=AMO+F%#( AP4—-AMO) 
CT=BMOtTF#(BP4~—-BNO0) 
RETURN 

Mert PLUSS (PD,AP8,8P8)} 
CALL PLUS (PD,AP4,BP4) 
F=(J~-.4)"%2.5 
CQOH=AP4SH+FH(APB-APG ) 
CT=BP4+F#(BPB-BP4G) 
RETURN 

END 
SeepeOUTINE PLUS] (PD,AP1,8P1) 
Meer. GE. »«63)GOFOQ10 
B= ~.4//*#PD-.015 
GOTOLL 
AP1=.00516*PD-.049 
BP]=-.4 

RETURN 

END 

SeORUUTINE PLUS& (PD,AP8,BPB) 
meee D.GE. .35)G0OTO10 
APB=—.06#PD 

GOTO11 

Mee. O2629D-. 03 

mee Del Teme4)COTO1L2 
meee o.GI. -4) GOTO1L3 
Peo =. 085#PD-.275 

RETURN 

BPB=.165#°D-. 243 

RETURN 
BP8=.243#PD-. 3382 
RETURN 

END 
SUBROUTINE PLUS4 (PD,AP4,BP4) 
PE(PD.GT. w4)GOTOLO 
PeCPD.LT. ~-%)GOTOL!L 
AP4=.015 

Gayoi2 
AP4=.0848#PD-.019 
GOTO12 
meee. | 25k OOo) 

Pet PD. L Vauee weno) Gls 
TRiPD 2Gl. S4cotTol4 
BP4=.245#PD-.9 
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RETURN 

BP4=.37#PD-.04 

RETURN 
peg—-.4)5#PD—. 158 
RETURN 

END 

SUBROUTINE ZERO 4 PD ,AMO,s BHO) 
Mme PD.b).-.4)GU0T010 
PrF(PD.GT. .4)GUTOLL 
AHMO=.0125 

GOTO12 
PMO=-.125#PD-.0375 
GOTO12 

peeo=. LOIVZ#PD-.0311 
Meee O.L!.-.625)C0TO1 3 
Mee PD.GT. 6-425) G0TOM 
BMO=.3108PD+.020 © 

me TtUR' 

me =.53ePD0+.1575 
RETURA 

P=. 53#PD-.0732 
RETURN 

FND 

PoeROUTINE MIN2 (PD,4M2,BM2) 
Merl Gf. eC45)G0TOLO 
MmeePD) .LT. —-0%75)G0TO11 
AM2=.9125 

GOTO12 

Pee =. 1 23*PD-.043 
GOTO12 

Mee —-—.132*PD-—.05 
eet0i2 

MmerrD LT. ~-625)G0TOL3 
feet) GT. ~-425)G0TO14 
2 —e310*#PD+.02 

RETURN 

Beg —=.53*PD+.1575 
RETURN 

Me = 2 53%PD—-.0/732 
RETURN 

END 

SUBROUTINE MING (PD,AM4,BM4) 
merPD.LT.~-425)G0I010 
Met PD .GTe. «275)60TO11 
AM4=.Q006 

GOTH12 

Pers =—.0925*PD-.0335 
Seroi2 
AM4=.9725*PD-.014 

Pe( PD) .Gl. eco ars 
PECPO «Ll. =. er eaters 
BM4=.184*PD+.076 
RETURN 

BM4=.5#PD-.202 

RETURN 

BM4=.4%#PD+.145 
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RETUR' 

END 

SUBROUTINE MIN8& (PD,AM8,BMB) 
EGP D wb bee. 55) head ome 
AMB=.07#PD 

GOTQ13 
AMB=-.028#PD-.054 
MmePD .LT. =345)6010M 
mae) «Gls <42) GOTNI? 
BM8=.09#PD+.2760 
RETURN 

BM8=.21#PD+. 33 


“RETURN 


BM8=.283#PD+.1652 
RETURN 

END 

SUBROUTINE MIQNE (PD,AML,BML) 
AM1=.064#PD-.018 

MeereD .GT. .6)60TO1} 
MerPD .LE. .400)G0TO10 
BML=-.88#?D+.743 
RETURY 

BM1=.391 

RETURN 

Pie =.215 

RETURN 

END 
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PUBROUTINE SHAPT Shiota) 
THIS SUBROUTINE COMPUTES SS orl On 


ENTERING ARGUMENT IS SHAFT ROTATION RATE (SN) IN REV PER SECOND 
PeePUl 1S FRICTION, 1ORRUE Ss Oi) eu; 


SN=ABS(SN) 
MetoN.lLE.~.417)GOT010 
QF=SN*6000. 

RETURN 

QF=2500. 

BerURN 

END 





APPENDIX & 


Ship Dynamics 


le Ship Propulsion Equations 
eae eee 
The ship propulsion equations with no external forces or 
. rudder movement are the two first-order, non-linear, differential 


equations given below: 


ay 


gee = = (thrust equation) 


dn 
it pepe ee Q, (torque equation) 
v = ship velocity (ft/sec) 
T = thrust (lbs) 
R = resistance (lbs) 
Lb.sec” 


m = mass of ship plus 10% —— ) 


rotational inertia of drive shaft Ibeft-sec“) 


i 


it 
Q = prine mover drive torque on propeller (lb-ft) 
Q. = shaft friction toraue (lb-ft) 
Q, = propeller torque (1b-ft) 
n = propeller angular speed (rev/sec) 
The two equations are coupled through the propsller terms T and 
Qn The equations were solved aivultancously on an IBM 1130 
computer using & fourth order Rungé-eKutta integration techniques. 
Ze Sign Convention 
The sign convention used here is that propelier angular spced, 
n, and velocity, v, are defined as positive for forward ship motion. 


Torques are considered positive when acting in the direction of 
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APPENDIX E (Cont'd) 


(Cont'd) 


positive ns; thrusts are considered positive in the direction of 


positive v. 


; Ship Displacement 


For al11 simulations, a ship with a displacement of 4000 tons 


was used. In the equations of motion, 10% of this figure was 


included to account for entrained water. 
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APPENDIX F 


FTYA~2 Drive Train 


Digital Simulation of an FT4A-2 Forward Drive Train 

Figure (Fel) is a sketch of a drive train powered by two 
FTtA-2 marine gas turbines that was used in the simulations. 
Except for the propeller the drive train is basically the same 
used by Rubis in Reference (11). The reduction gear ratios relate 
propeller speed and torque to power turbine speed and torque. 

The drive train moment of inertia was assumed to be the same 
as Rubis', Wize, AGO x 10° lb-ft~sec“ (referred to propeller 
shaft). 

it should be noted that at high speeds and low fuel rates, the 
FT4A«2 gas turbine subroutine will yield negative toraues. A 
statement was made in the main program to ignore this and sub- 
stitute instead a linear windage torque of the form, 


=? + 3000 (ft-lb) 


ti 


Q 

N = Free Turbine Speed (RP#H) 

For the crashback simulation, a Falk clutch whose torque 
capacity is 

Q, = 5687.5  Pnet 
where 

Phet = net air pressure acting on shoes 
was used (cf. section 3.2). A supply air pressure of 200 psi was 
assumed. The criteria for clutch disengagement was the same used 
by Rubis in Reference (11): whenever the propeller toraue "seen" 


by the clutch exceeded clutch torque, clutch disengagement was 
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APPENDIX F (Cont'd) 


(Cont'd) 
was assumed to have taken place instantaneously. For the crashe- 
back simulation, the equation for forward clutch disengagement 


torque was 


Q, = 568765 (195-30tuN*) (£t-2b) 
t = time (sec) 
N = propeller shaft rotation rate (rev/sec) 


The propeller torque opposing the forward clutch was 


. a s 
po _pP... = 
“pe © =k G35197 


3 


Ks = reduction gear ratio (see fig. (F-1)) 


°p 


During reverse clutch engagement, reverse clutch torque is given 


i 


propeller torque ft-lbs 


by 
Qop = 5687.5 (St~S.lin”) (£t-1b) 
1 = time (sec) 
n = propeller shaft rotation rate (rev/sec) 


The open drive train is defined to exist whenever the auantity 
(5t-Se4n@) is negative; if this occurs, Que is set to a zero 
value obviating a negative value for clutch torque. <A slip speed 
across the clutch exists during reverse clutch engagement and 


this is given by 


Nelip = K 3 ~ (= é 
Ce 


propeller rotation rate (rpm) 


Lamm 
a4 
ih 


po 
it 


3 power turbine speed (rpm) 
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APPENDIX F (Cont'd) 
(Cont'd) 


Ky s Kos Kas = reduction gear ratios 


int ° a ; a i . see | 
Tne product of Nsiip and QR is the power dissipated in the clutch, 


' and the integral of power over time is a measure of the energy 


absorbed in the clutch. 
Reverse clutch engagement was assumed to exist the instant 


W equaled zero. 


slip 
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APPENDIX G 


IM2 500-A Drive Train 


Digital Simulation of an LiM2500-4 Powered Drive Train 


Figure (G-1) is a sketch of the drive train powered by two 
Li2500eA marine gas turbines with a controllable reversible pitch 
propeller (CRPP). A popular control strategy is to increase speed 
by first increasing propeller pitch while keeping shaft rotation 
rate constant; once full pitch is attained, shaft rotation rate 


is increased.(24) A proportionaleintegral controller of the form 


Wp = 54 235 
RPI 


ii 


We 


RPM = (ordered shaft RPM)-(Actual Shaft RPM) (RPM) 


fuel rate (lb/hr) 


S 


differential operator 


ft 


was used in one simulation to keep shaft rotation rate constant 
while increasing pitch. 

As in the FTA-~2 subroutine, the LM2500-A subroutine will 
also yield negative torques for high power turbine RPM and low 
fuel rates. These torques were also ignored end a windage torque 
expressed by 


Cy 


N 


-1.965 N + 5500 (ft~-1b) 


i 


power turbine speed (RPM) 
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SUBROUTINE EQS RE 

ACCELERATING FRGM. SRKNG TS Syl Tice 1 4teeZ 

COMMON Tye OTs YU2ZO) se 4 200% SIME ee ee 
VA=.96*Y(2) 

CALL PRPLR (VA,CT;CO) 

DENOM=VA##24+(33.8Y(1) )##2 

COMPUTE PROPELLOR TERGBEyY(5} 
Y¥(5)=262.5*#DENOM#CG 

eeMPUTE PROPELLOR THRUST 

Y(6)=171.5#DENODM&4CT 


Meee OPELLOR THRUST IS REDNCED BY THRUST DEDUCTION FACTOR 


memo} =Y(6)#,98 

foe 1S REDUCTION GEAR EFFICIENCY AS A FUNCTIGN OF SitART S8ErD 
Me. (9+. 22%(1.0-EXP4-=Y(1)/71.04)) 

fame 1S FUEL FLOW AYA FUNCTION OF TIME FOR THIS =VANEU Ye 
ee) =2900.+(1.-EXP(—71 )) #9500, 

Peboe oUIP RESTS VANCES oe Ve EUNCTIGN OF “SHIP *S SPEER 
RT=104.8#(Y(2) ) #2 

Pee | COMPUTES SHAFT FRICTION 

SN=Y(1) 

eet SHAFT (SN,CF) 

PNS=861.%#Y(1) 

WF=Y(7) 

pepe COMPUTES POWER TURBINE TORQUE 

CALL FT4A2 (AN3,WF e003) 

MA Y=.3+.7#(1.-EXP(-.45T) ) 

Q3=03-814. 

QO3=OQ38DEL AY 

Q3=03+814. 

¥«9) IS POWER TURBINE TCRQUE 

Y¥(9)=93 

QD IS POWER TURBINE TORQUE MODIFIED BY REDUCTION GEAR FACTOR 
QD=28. 728 Q3 ETA 

femme lS SHAFT SPEED IN RPS 

fee = ( QOD-QF-Y¥(5))#5.54E-7 

tome) IS SHAFT SPEED IN RPM DIVIDED BY 100. 

Y¥(8)=.6#Y(1) 


fom iS SHIP'S VELOCITY IN FPS 


Beeice )=(Y(G)-RT)43.33E-6 

fee 1S DISTANCE TRAVELLED BY SHIP 
met 3)=Y(2) 

fee? 1S SHIP'S VELOCITY IN KNOTS 

MAS )=.591#Y(2)} 

fomro) 1S PROPELLER THRUST COEFFICIENT 


flo) =CT 

feet) [S PROPERDESR TORQUE COEFFICIENT 
Y(17)=CG 

mea?) «GT. LOsIDV=1., 

RETURN 


END 





DAGE 


) 


140 


l 05/7032 


SUBROUTINE EQSIM 

STOPPING THE SHIP WITH THE FTI4A-2 

COMMON Ty DOT,Y(20),DY(20) »,STIME,FTIME»NEWDT 
VA=.90#Y( 2) 

GALL PRPiLR Ae eG een) 
DENGM=VA##24(33.#Y(1))##2 

COMPUTE PROPELLOR TORQUE Y(5) 
¥(5)=262.5*#DENOMN#CQ 

COMPUTE PROPELLOR THRUST 
Y(6)=1L71.5#DENOMN#CT 

PROPELLOR THRUST JS REDUCED BY THRUST DEDUCTION FACTOR 
Y(6)=Y(6)#.98 

ETA IS REDUCTION GEAR EFFICIENCY AS A FUNTION OF SHAFT SPEED 
ETA=.75+.22#(1.0-EXP(-ABS(Y(1))/1.04)) 

Y(7) IS FUEL FLOW AS A FUNCTION OF TIME FOR THIS MANFUVER 
reer = 1600.4+2600.#(T-115119) 

meee) «Gl. 12000.)Y¥{7T)=Si2000. 

RT IS SHIP RESISTANCE AS A FUNTION QF SHIP'S SPEED 
RT=104.8#(Y(2))##2 

SHAFT COMPUTES SHAFT FRICTION 

= Y (1) 

CALL SHAFT (SN,QF) 

AN3=600.25#(-Y(1)) 

Y(9) IS POWER TURBINE SPEFD IN RPM 

Y(9)=AN3 

WFE=Y(7) 

FTI4A2 COMPUTES GAS TURBINE TORQUE 

CBLL FI4A2 (AN3,WF,O3) 

MrP1Q3 elles 0. )Q3=-2.4AN34+3000. 

Q3=03-2165.2 
DELAY=.3+.7#(1.-EXP(-.4#(T-11.119))) 
Q3=Q3#DELAY 

@2=03+2165.2 

Y(8) IS POWER TURBINE TORQUE DIVIDED BY 100 
¥(8)=Q3#.01 

QD IS GAS TURBINE TORQUE MODIFIED BY REDUCTION GEAR FACTOR 
QD=-20.1*Q3#ETA 

Y(1) IS SHAFT SPFED IN RPS 
DY(1)=(QD-QF-Y(5))#6.75E-7 

ce) 6ULS SIME * SeVEMOCITY IN FPS 
DY(2)=(Y(6)-RT)#3.33E-6 | 

Y(3) IS DISTANCE TRAVELLED BY SHIP 

mrt 3)=Y(2) 

m4) 1S SHIP SS VElOCITY IN KNOTS 

re 4)=.591#Y( 2) 

Y(19) IS THRUST COEFFICIENT 

m1 9) =CT 

¥(20)= IS TORGUESERIEFPFICIENT 

yeZ0)=CC 

PFCY(2).CEf Cee =i 

Li {1 .Gl. owen 

RETURN 

END 
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APPENDIX H 


KV Major 12 Drive Train 


Figure (H-1) is a sketch of the drive train used for simulation. 


The WR° of the engine-flywheel was provided by the Hawker-Siddeley eles 
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APPENDIX I 


B&W 7K9CFF Drive Train 


Figure (I-1) is a sketch of the 7K98FF drive train. The WR° 


of the engine was not available in the Dieemetrnae surveyed. An 
official of the Burmeister and Wain Corporation previded the following 
weights s 

Piston eccee 418 tons 

Crossn@adece 3635 tons 

Crank eccce 4097 tons 
From the Shock and Vibration Handbook", Vol. 3, po 38~2, the following 


expression was found for finding the wR of a diesel: 


y= (op + W, (le dR 

me = weight of piston, piston pin, and cooling fluid (ib) 

W, = weight of connecting rod (1b) 

h = fraction of rod length from the lower center of gravity 
R = crank radius (in) 


From the above expression, the polar moment of inertia was estimated 


mee X 10? Lb-ftesec“, 
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